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ABSTRACT 

We derive lower bounds on intergalactic magnetic fields (IGMFs) from upper 
limits on the pair echo emission from the blazar Mrk 501, that is, delayed GeV 
emission from secondary e~e + pairs produced via interactions of primary TeV 
gamma rays with the cosmic infrared background. Utilizing only simultaneous 
GeV- TeV light curves observed by VERITAS, MAGIC and Fermi-LAT during a 
multiwavelength campaign in 2009 that included a TeV flare, bounds are deduced 
on the IGMF strength of B > 10~ 20 G at 90% confidence level for a field coherence 
length of 1 kpc. Since our analysis is based firmly on the observational data alone 
and nearly free of assumptions concerning the primary TeV flux in unobserved 
periods or spectral bands, our evaluation of the pair echo flux is conservative and 
the evidence for a non-zero IGMF is more robust compared to previous studies. 

Subject headings: BL Lacertae objects: individual (Mrk 501) — gamma rays: 
observations — gamma rays: theory — magnetic fields — radiation mechanisms: 
nonthermal 

1. Introduction 

Although their existence is yet to be observationally confirmed, the possibility that 
weak magnetic fields were generated ubiquitously in the early Universe has attracted 
considerable attention. Such "cosmological" magnetic fields are potentially important from 
at least two perspectives. First, by serving as the seed fields for subsequent amplification by 
galactic dynamo mechanisms, they may have been the ultimate origin of the magnetic fields 
seen today in galaxies and clusters of galaxies (Widrow 2002). Second, in some regions 
such as the centers of intergalactic voids, they may have survived to the present day as 
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intergalactic magnetic fields (IGMFs) without being affected by later magnetization from 
astrophysical sources (Bertone et al. 2006), and therefore may provide us with valuable, 
fossil information about physical processes in the early Universe. So far, various mechanisms 
have been proposed for the generation of such cosmological magnetic fields, with predicted 
field amplitudes in the range 10~ 25 - 10~ 15 G (e.g. Gnedin et al. 2000; Langer et al. 2005; 
Takahashi et al. 2005; Ichiki et al. 2006; Ando et al. 2010). Although possibly sufficient 
as seeds for galactic dynamos, such tiny magnetic fields are extremely difficult to confirm 
observationally through conventional methods such as Faraday rotation measurements or 
their effect on cosmic microwave background anisotropics (Widrow 2002). 

In this context, a potentially powerful probe of weak IGMFs with strengths 
10~ 20 — 10~ 15 G may be offered by pair echo emission from extragalactic TeV sources such as 
blazars or gamma-ray bursts (e.g. Plaga 1995; Dai et al. 2002; Razzaque et al. 2004; Ichiki 
et al. 2008; Murase et al. 2008; Takahashi et al. 2008, 2011; Neronov & Semikoz 2009). 
Pair echos comprise inverse Compton (IC) emission from secondary e"e + pairs produced 
via intergalactic 77 interactions among primary TeV gamma-rays and infrared-UV photons 
of the extragalactic background light. They are distinguishable through their characteristic 
time delay and spectral variation that depend on the properties of the intervening IGMFs, 
while being insensitive to galactic-scale magnetic fields either local to the source or the 
observer. An alternative approach focusing on the spatial extension of the secondary 
gamma-ray halo emission has also been discussed (Neronov & Semikoz 2009; Ando & 
Kusenko 2010). 

Recently, Neronov & Vovk (2010) claimed a lower bound on the IGMF of order 10~ 15 G 
from the non-detection of pair-echo components in the GeV spectra of selected TeV blazars. 
A few other studies also gave similar results (Dolag et al. 2011; Tavecchio et al. 2010, 
2011). However, as pointed out by Dermer et al. (2011), an implicit but crucial assumption 
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in deriving these bounds was that the TeV emission has been persistent for at least the past 
10 6 years, at the level observed in these objects on a small number of specific dates in the 
last several years. This assumption is quite questionable, because TeV blazars are generally 
known to be highly variable, with their TeV flux fluctuating by more than two orders of 
magnitude over timescales of several years and less, for sufficiently well observed objects 
(e.g. Bartoli et al. 2011 and references therein). With a more relaxed assumption that TeV 
gamma-rays from the blazar 1ES 0229+200 was steady only during the years 2005 to 2009 
when the observations were conducted, Dermer et al. (2011) obtained a much weaker lower 
bound of order 10~ 18 G. In fact, even this assumption is not entirely satisfactory, since the 
actual time coverage of TeV observations in this period was very sparse, being only several 
days in 2005, 2006 and 2009. Taylor et al. (2011) have presented a similar analysis for a 
few blazars with analogous assumptions on their TeV activity during unobserved periods. 

Here we constrain the pair echo component from Mrk 501 and derive lower bounds 
on the IGMF, relying solely on observational data from a multiwavelength campaign over 
a few months in 2009. Thanks to a flare detected by VERITAS during this period, lower 
bounds of order 10~ 20 G can be derived, as the expected pair echo flux for weaker IGMFs 
would exceed the concurrent upper limits on the daily GeV flux from Fermi-LAT. Since we 
do not impose any assumptions on the TeV emission preceding the campaign, our bounds 
implying non-zero IGMFs are much more robust than those obtained previously. 

2. Mrk 501 observations 

Abdo et al. (2011) reported on a 4.5 month-long multiwavelength observing campaign 
of Mrk 501 in 2009, featuring TeV data from VERITAS and MAGIC, and GeV data 
from Fermi-LAT. The TeV spectra and light curves are essential input for evaluating the 
expected properties of the pair echo, while the GeV data constrains the echo itself that 
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typically emerges in this energy range (§3). 

The VERITAS observing runs were conducted intermittently over a three month 
period MJD 54907-55004, resulting in a total of 9.7 hours of good quality data. During 
the 3-day "flare state" of MJD 54953-54955, the TeV flux increased significantly, by about 
a factor of five relative to the remaining period, referred to as the "quiescent state". 
The measured spectra for either state can be fitted by power-law functions of the form 
log F(E) = log isT - alog(£/TeV), with K = (4.17 ± 0.24) x 10- n ph/cm 2 /sec/TeV and 
a = 2.26 ± 0.06 for the flare state, and K = (0.88 ± 0.06) x 10- n ph/cm 2 /sec/TeV and 
a = 2.48 ± 0.07 for the quiescent state. MAGIC was not available during the flare state but 
gave a spectrum consistent with VERITAS for the quiescent state. Both states are taken 
into account for calculating the pair echo emission. 

On the other hand, Fermi-LAT performed uninterrupted monitoring of Mrk 501 in the 
survey mode during MJD 54683-55162, which included the whole campaign. Although the 
GeV data did not reveal any strong flaring activity, some mild flux variations were detected 
on timescales of about 30 days. In contrast, changes in the contemporaneous GeV spectra 
were conspicuous, particularly the 30-day spectrum for MJD 54952-54982 that was much 
harder than other periods. Note that the first three days of this span correspond to the TeV 
flare state, which may have possibly lasted much longer but was missed by the sparse and 
irregular time coverage of TeV observations. Nevertheless, below we conservatively assume 
the flare duration to be just these three days. 

The Fermi-LAT data was acquired via the Fermi Science Support Center (FSSC) |*j 
and analyzed to search for day-timescale flux variations around and after the flare. We 
use the standard analysis tools supplied by FSSC in three energy bands, 100 MeV-1 GeV, 
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1 GeV-10 GeV and > 10 GeV, so as to keep reasonable photon statistics in time intervals 
as short as one day while retaining some energy resolution. Since the statistics is still small, 
we adopted the aperture photometry method, i.e., events falling within one degree from the 
source were counted. Note that background events above 1 GeV in one-day bins at the high 
Galactic latitude of Mrk 501 are essentially negligible. Flux upper bounds were calculated 
assuming Poisson statistics when the data show no significant gamma-ray signal. 



3. Pair echos and application to Mrk 501 

First we briefly summarize the basic physics of pair echo emission (see e.g. Ichiki et 
al. 2008 for more details), and also provide an improved formulation for application to 
Mrk 501. Primary gamma-rays with energy E 1 > 1 TeV emitted from an extragalactic 
source have mean free path A 77 = 1/(0.26<tt^ir) = 190 Mpc (nrR,/0.01 cm -3 ) -1 for 77 pair 
production interactions with photons of the cosmic infrared background (CIB), where a T 
is the Thomson cross section and n IR is the number density of CIB photons most relevant 
for the interactions. The produced pairs with energy E e E 1 /2 give rise to the pair echo 
emission by IC upscattering of ambient cosmic microwave background (CMB) photons 
to average energy (E echo ) = 2.7T CMB 7e = 2 - 5 GeV ( E i/ 2 Te V) 2 , where 7 e = E e /m e c 2 
is the Lorentz factor of the pairs and T C mb = 2.7 K is the CMB temperature. Thus, 
primary gamma-rays in the range E 1 ~ 1 — 5 TeV induce echos with typical energies 
E echo ~ 1 — 10 GeV. The IC mean free path of the pairs is Aic, SC at = 1/(°"t^cmb) = 1-2 kpc, 
where wcmb ~ 420 cm~ 3 is the CMB photon number density. The pairs upscatter CMB 
photons successively until they lose most of their energy after propagating an IC cooling 
length A IC ,cooi = Zm 2 J(AE e a T U CMB ) = 350 kpc (E e /1 TeV)" 1 , where U CMB is the CMB 
energy density. The length scales for A 77 and Aic,cooi imply that the secondary pairs 
typically arise in locations far removed from the source on scales of intergalactic voids, 
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whereas the pairs propagate only for short distances within such regions while generating 
the echo emission. 

The pair echo emission arrives at the observer with a time delay rela- 
tive to the primary emission, caused by the effects of angular spreading in pair 
production and IC interactions, as well as by deflections of the pairs in in- 
tervening magnetic fields. The typical delay time due to angular spreading is 
At ang = (A 77 + Aic,cooi)/27 e 2 ~ 3 x 10 3 sec (E ccho /l GeV)- 1 (n m /0.01 cm^t' (Ichiki et 
al. 2008), while that due to magnetic deflections is At B = (A 77 + AiccooiX^l)/^ where 
(f?!) 1 / 2 = max [Aic,cooi/ r L, (Aic,cooi r coh/6) 1/,2 /rL] is the variance of the magnetic deflection 
angle in terms of the Larmor radius r L and field coherence length r coh . If r coh <C Aic,cool> 
At B « 2 x 10 4 sec (E echo /l GeV)-V 2 (B /lO" 19 G) 2 (r coh /l kpc)(n IR /0.01 cm" 3 )- 1 , where 
B is the field amplitude. Hereafter we fiducially take r co h = 1 kpc (see e.g. Langer et al. 
2005), but the results are scalable to other values of r co h since it enters only through the 
combination B 2 r coh and only when r coh < Aic, C ooi- The total delay time is approximately 
At = At ang + AtB, and the magnetic field properties are reflected in the delay as long as 
At ang < At B . 

The spectra and light curves of the pair echo can be evaluated as follows. For a primary 
fluence dN 7 /dE y , the associated time-integrated flux of secondary pairs is 

^( 7e ) = 4m e ^I(£ 7 = 2m e7e ) [l - e -^= 2 >-)] , (1) 

where r 77 (i? 7 ) is the 77 optical depth in the CIB. The time-dependent echo spectrum is 

d 2 N echo f dN e d 2 N lc 

, TTTET , ( 2 ) 



dtdE 7 J c?7 e dtdE, 

where d 2 Nic/dtdE-y is the IC power from a single electron or positron, and dN e /d^ e is the 
total time-integrated flux of pairs responsible for the echo emission observed at time t, 
related nontrivially to dN e fl/d^ e in Eq. (pQ). If the distance to the source from the observer 
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D 3> A 77 , it can be approximated by dN e /di e = (Xic, C ooi/cAt)dN e;0 /di e (Dai et al. 2002), 
but this is not the case for Mrk 501 whose D ~ 130 Mpc can be comparable to A 77 for 
E 7 > TeV. In order to improve the evaluation of dN e /di e by accounting for the finite 
probability of pair production near the observer, we consider the time-integrated flux due 
to the fraction of pairs that originate between radii r 77 and r 77 + Ar from the source, 

= 4m e ^ ( e - r WA 77 _ e -(r 77 +Ar)/A 77 ) (3) 

The total time-integrated flux of pairs can be evaluated by integrating over r 77 as 

^ = r*„_^=u m .^fpi i ( 4 ) 

d^e J cAt(r 77 ) dE^ A 77 

where At(r 77 ) is given by the above expression for At = At ang + Ate with A 77 replaced by 

r 77 . 

Note that the pair echo fluence is determined by the total amount of absorbed primary 
gamma rays and thus independent of the IGMF, in contrast to the pair-echo flux which 
is roughly given by the fluence divided by At. Weaker IGMFs generally give higher echo 
fluxes, as long as the time delay does not become dominated by angular spreading and the 
echo flux remains sensitive to B. For r coh = 1 kpc, At B approaches At ang if B ~ 10~ 20 G. 

In applying the above formalism to the 2009 Mrk 501 activity, we clarify what we 
employ for the primary TeV spectra and light curves. Both flare and quiescent states can 
make important contributions to the pair-echo emission. A TeV flare was observed for at 
least three days from MJD 54953 (§2); however, it may have continued for a longer time, 
or even separate flares could have occurred over the following weeks, as can be speculated 
from the hard, 30-day Fermi-LAT spectrum. Nevertheless, we choose to be conservative 
and assume that there is no other flare state during the campaign besides the three days 
seen by VERITAS. Although the quiescent state was also only sparsely sampled at TeV, 
since both VERITAS and MAGIC measured a consistent flux and spectrum at separate 
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times, we postulate that the quiescent emission is steady over the period covered by the 
TeV telescopes, the sole assumption we make regarding TeV activity not directly observed. 
Thus, we consider the primary light curve to consist of a flare state with a top-hat shape 
for the 3 days MJD 54953-54955, together with a steady quiescent state for the preceding 
46 days MJD 54907-54952 as well as the ensuing 49 days MJD 54956-55004. The primary 
flux and spectrum for each state are chosen such that they are compatible with those 
observed by VERITAS after accounting for intergalactic 77 absorption with the CIB model 
of Franceschini et al. (2008), and are described by the same power-law functional form 
mentioned in §2 but with the parameters K = 9 x 10 _11 ph/cm 2 /sec/TeV and a = 2.0 for 
the flare state, and K = 2 x 10 _11 ph/cm 2 /sec/TeV and a = 2.3 for the quiescent state. 
Minimum and maximum cutoffs are also imposed at 0.1 TeV and 5 TeV, respectively, the 
latter corresponding to the highest energy photons detected by VERITAS and MAGIC. 
Comparing the pair echo emission calculated in this way with the observed GeV limits gives 
conservative lower bounds on the IGMF, since any additional primary emission, outside 
either the above time interval or the above spectral range, would only add to the pair echo 
flux and lead to tighter bounds. 

4. Results and discussion 

Fig. [1] shows the spectra of the primary and pair echo emission for the flare and 
quiescent states when B = 10~ 20 G. The primary spectra are displayed both with and 
without intergalactic 77 absorption, the latter to be compared with the absorption-corrected 
VERITAS data as given in Abdo et al. (2011). The echo from the flare state is plotted 
at observer times t =1, 10 and 100 days after the flare, fading progressively on timescales 
approximately corresponding to At. In contrast, here the echo due to the quiescent 
state is essentially stationary on the timescale of the campaign and independent of B. 
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Note, however, that for stronger B with accordingly longer Ate, even the quiescent echo 
component can be nonstationary, particularly at low energies. Only when the primary 
emission persists at a steady level for a time considerably longer than does the echo 
reach stationarity, as demonstrated by Dermer et al (2011). 
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Fig. 1. — Spectra of the primary and pair echo emission of Mrk 501 for the case B = 10~ 20 G. 
The primary spectra for the flare and quiescent states are each shown with (long-dashed) 
and without (solid) intergalactic 77 absorption, along with the absorption-corrected data 
from VERITAS observations. Also plotted are the echo from the flare state at observer time 
t =1, 10 and 100 days after the flare (dashed, from top to bottom), as well as the echo from 
the quiescent state (dotted). 



The light curves of the pair echo in the 1-10 GeV band after the onset of the TeV flare 
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on MJD 54953 are plotted in Fig. [2] for different values of B. The flux initially rises over 
the duration of the primary flare, and then decays roughly exponentially on timescales of 
order At as the primary emission switches to the quiescent state. At sufficiently late times, 
only the quiescent emission contributes to the echo and its flux approaches a steady level. 
For weaker B, the echo flux responds to changes of the primary flux more quickly by virtue 
of the shorter Aig, and the peak flux is larger and more susceptible to observational limits. 

Thus, in order to observationally constrain IGMFs in the range B ~ 10~ 20 — 10~ 19 G, 
GeV-band light curves with time resolution of order a day are necessary. We analyzed 
the Fermi-LAT data to obtain 1-10 GeV gamma-ray fluxes and upper limits with 1-day 
time binning, the results of which are compared with the pair echo expectations in Fig. [2j 
Most are upper limits, which is not surprising given the limited sensitivity of Fermi-LAT 
with such short integration times. The strongest limits come from the second, third and 
fourth days, being comparable to the pair-echo predictions for B < 10~ 20 G and providing 
important lower bounds on IGMFs in this range, while higher values of B cannot be usefully 
constrained by the current analysis. Note that primary emission in the GeV range is also 
generally expected and should contribute significantly to the LAT light curve, which implies 
that the actual pair echo flux is even lower and the true IGMF lower bounds stronger. 
However, given the lack of reliable knowledge on the primary GeV spectra and variability, 
we restrict ourselves to conservative constraints by not accounting for any such components. 

To deduce bounds on B, we first calculate from the obtained LAT upper limits the 
probability for each day that the predicted pair-echo flux for a specific value of B does not 
exceed the true flux. Then we regard each value of B as being allowed at the probability 
equal to the product of the above daily probabilities. Thus we arrive at our main result 
that B > 10 -20 G at about 90% confidence level, determined mostly by the limits from 
the second through fourth days. We have also carried out similar analyses for other energy 
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Fig. 2. — Light curve of Mrk 501 in the 1-10 GeV band from the onset of the TeV flare on 
MJD 54953. Pair echo expectations for B = lO" 205 , lO" 20 , 10~ 19 - 5 , 10~ 19 G (curves from top 
to bottom) are compared with Fermi-LAT data binned at 1-day intervals using the aperture 
photometry method, where errors (vertical bars) or upper limits (downward arrows) are at 
68% confidence level. 

bands < 1 GeV or > 10 GeV but were not able to derive significant constraints, as can be 
expected from Fig. CD 

In summary, by comparing the daily GeV flux upper limits from Fermi-LAT for the 
blazar Mrk 501 during and after its TeV flare in 2009 with the expected light curves for 
the the associated pair echo emission, we have derived lower bounds on IGMF strengths of 



- 13 - 



B > 10~ 20 G at 90% confidence level for a field coherence length r coh = 1 kpc. The result 
can be roughly scaled for other values of r coh as B > 5 x 10~ 22 max[(r coh /350 kpc)" 1 / 2 , 1] G. 
This bound is weaker compared to other recent results obtained through similar methods, 
which, however, all relied on unproven assumptions regarding the TeV emission during 
unobserved periods on timescales of years (Dermer et al. 2011; Taylor et al. 2011) or much 
longer (Neronov k Vovk 2010; Dolag et al. 2011; Tavecchio et al. 2010, 2011). In constrast, 
our analysis is entirely free of such assumptions other than for the quiescent state, and thus 
can be considered the most robust indication so far for the existence of non-zero IGMFs. 

Future progress can be expected from regular, long-term coverage of the multi-TeV 
emission by wide-field facilities such as HAWC or LHAASO, as well as GeV-TeV spectral 
variability measurements with high sensitivity and time resolution by CTA to disentangle 
and positively identify the echo component from the primary emission. If IGMFs are 
somewhat stronger, the spatially-extended pair halo emission may be detectable by CTA. 
Such advances will surely pave a new road toward understanding cosmic magnetic fields. 

We thank Andrii Neronov for crucial comments. This work is supported by Grants- 
in-Aid from MEXT of Japan, No. 21840028 (KT), No. 22540315 (MM), No. 21740177, 
No. 22012004 (KI) and No. 22540278 (SI), and for the global COE program "The Next 
Generation of Physics" at Kyoto University and "Quest for Fundamental Principles in the 
Universe" at Nagoya University. 
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